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Разработан алгоритм, имитирующий процесс зарождения и 
развития усталостных трещин в сварном соединении на примере 
узла крепления мобильной платформы и опорной обечайки 
поворотного устройства грузоподъемного крана бой\уа НМК 
170 ЕС (63 т), разрушение которого произошло в 2011 году на 
территории ОАО «Новороссийский морской торговый порт» 
(НМТП). Предложенный алгоритм учитывает неоднородную 
геометрию шва, влияющую на значение коэффициента 
концентрации напряжений, и неравномерное распределение 
номинальных напряжений, связанное с конструктивными 
особенностями узла. В основе процесса моделирования лежит 
экспериментально установленный факт, что профиль геометрии 
перехода от шва к основному металлу вдоль сварного 
соединения изменяется случайно. В результате этого появление 
трещин вдоль шва носит случайный характер. При циклическом 
нагружении развитие разрушения представляет собой 
последовательный процесс зарождения трещин на различных 
участках сварного соединения вдоль шва, подрастание 
одиночных поверхностных трещин по законам механики 
разрушения, появление новых трещин, слияние трещин на 
соседних участках. При слиянии берегов трещины изменяется 
профиль контура, что приводит к увеличению скорости роста 
глубины трещины. В дальнейшем процесс развития 
многоочагового разрушения продолжается до достижения 
трещиной критического размера. Скорость накопления 
повреждений и развития разрушения зависит от степени 
неоднородности концентрации напряжений вдоль шва. 
Реализация предложенной модели позволит решить задачи 
определения периодичности диагностирования узла, мест с 
наибольшей вероятностью возникновения усталостных трещин 
и обосновать возможные методики усиления конструкции. 


Ключевые слова: многоочаговое разрушение, концентрация 
напряжений, распределение НДС, неравномерное 
распределение номинальных напряжений. 


Введение. Для прогнозирования процесса разрушения сварных конструкций при циклическом нагружении зачастую 
используют аналитические модели, которые позволяют оценивать долговечность и надежность относительно простых 
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тел с трещинами. Однако они не всегда позволяют учитывать многообразие факторов, влияющих на процесс 
разрушения сварных соединений, и вынуждают прибегать к упрощению реальных явлений. Этот факт указывает на 
то, что использование аналитических моделей для расчета характеристик прочности сварных соединений может быть 
слишком неточным. 

Имитационное моделирование позволяет учитывать большое количество особенностей процесса развития 
разрушений сварных соединений, особенно в крупногабаритных конструкциях, которые могут иметь несколько зон 
концентрации напряжений. 

Процесс статистического моделирования зарождения и развития разрушений в сварных конструкциях 
подробно описан в работах сотрудников кафедры «Машины и автоматизация сварочного производства» ДГТУ [1-3]. 
Однако в этих работах моделирование рассматривалось с позиции неизменного распределения номинальных 
напряжений вдоль шва. 

Целью работы является обоснование возможности использования процесса имитационного моделирования с 
учетом неоднородного распределения номинальных напряжений в конструкции. 

Как показали результаты многочисленных экспериментов [1-3], зона разрушения конструкций с 
протяженными угловыми швами связана с концентрацией напряжений в месте перехода от шва к основному металлу. 
При этом кинетика развития разрушения существенно зависит от неоднородности геометрии профиля шва вдоль 
протяженных сварных соединений. Для сварных соединений с однородной геометрией профиля шва вдоль 
соединения разрушение происходит в виде появления одиночных, далеко расположенных друг от друга трещин и их 
роста в соответствии с законами механики разрушения для одиночных поверхностных трещин. В этом случае 
долговечность конструкции с учетом стадии развития трещины можно рассчитать по следующей формуле: 

М=М,+М,, (0 


где №, — длительность периода зарождения трещины; № „—_ длительность периода роста одиночной трещины [4]. 


С увеличением степени неоднородности геометрии профиля шва вдоль сварного соединения количество 
центров зарождения трещин возрастает. Фронты отдельных трещин могут соприкасаться, что приводит к слиянию 
трещин в процессе их роста при циклическом нагружении. Такое разрушение получило название «многоочаговое». 

В отличие от случая одиночных или редко расположенных трещин, кинетика процесса многоочагового 
развития представляет собой параллельное развитие трех процессов [3]: 

— зарождение одиночных трещин на различных участках сварного соединения; 

— рост одиночных трещин до соприкосновения берегов трещины с трещиной на соседнем участке; 

— слияние соседних трещин с образованием единого фронта трещины большей протяженности. 

При дальнейшем циклическом нагружении процесс роста существующих трещин, образования новых трещин 
и их слияние продолжается до достижения одной из трещин критического значения. 

Для случая многоочагового развития разрушения использовать формулу 1 не представляется возможным, т. к. 
при слиянии берегов соседних трещин форма контура объединенной трещины существенно изменяется и для 
моделирования кинетики такого развития разрушения необходимо разрабатывать имитационные модели. 

В работе [5] доказано, что расположение участков сварного соединения, на которых зарождаются трещины, 
должно носить случайный характер, поскольку зарождение трещин зависит от амплитуды локальных деформаций. 

В соответствии с уравнением Басквина-Мэнсона-Коффина [6-8] длительность периода зарождения 
разрушения при циклическом нагружении определяется амплитудой деформаций в локальной зоне сварного 
соединения. Величина локальной амплитуды деформаций зависит от уровня номинальных напряжений и 
коэффициента концентрации напряжений в зоне перехода от шва к основному металлу. 

Номинальные напряжения в зоне расположения сварного шва зависят от конструкции рамы платформы, 
опорно-поворотного узла и нагрузки на стрелу крана. В процессе поворота стрелы крана уровень номинальных 
напряжений в зоне сварного соединения изменяется, что создает неравномерное распределение их вдоль шва. 

Для анализа напряженно-деформированного состояния (НДС) металла платформы в зоне сварного 
соединения при различных углах поворота стрелы использовался метод конечных элементов (программный продукт 
АМЗУ5). 

Результаты расчета представлены на рис. 1, где пунктирной линией показано распределение номинальных 
напряжений вдоль сварного соединения при расположении стрелы вдоль оси Г-Г; штрихпунктирной линией — при 
расположении стрелы вдоль оси П-П; сплошной линией — при расположении стрелы вдоль оси Ш-Ш. Независимо от 
угла поворота стрелы номинальные напряжения достигают максимальных значений в одних и тех же точках (зона 1; 
зона 2; зона 3). 

Анализ конструкции рамы показывает, что эти зоны соответствуют местам расположения ребер жесткости 


2 рамы, которые изменяют локальную жесткость конструкции [9]. 
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Рис. 1. Эпюра номинальных напряжений вдоль сварного шва, соединяющего раму 
платформы и обечайку опорного кольца 


Таким образом, в отличие от известных ранее моделей [3], в данном случае при разработке имитационной 
модели необходимо учитывать неравномерный характер распределения номинальных напряжений вдоль шва. Кроме 
того, проведенный конечно-элементный анализ показывает, что путем изменения конструкции рамы можно 
существенно влиять на ресурс работы сварного соединения [10]. 

Следующим фактором, определяющим амплитуду деформаций при циклическом нагружении сварного 
соединения, является коэффициент концентрации напряжений в зоне перехода от шва к основному металлу. 

Основные параметры, оказывающие влияние на коэффициент концентрации напряжений тавровых сварных 
соединений, проиллюстрированы на рис. 2. 








Рис. 2. Геометрические параметры таврового соединения 


В работе А. В. Ильина, Г. П. Карзова и В. П. Леонова [11] для определения коэффициента концентрации 
напряжений тавровых сварных соединений использовалось полученное авторами опытным путем соотношение: 
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Наибольшее влияние на коэффициент концентрации напряжений (4. оказывает соотношение 7/5. Так как $ 





является постоянной величиной, то для определения а, главным фактором является определение радиуса перехода г. 
Он не регламентируется нормативными документами и может варьироваться от сотых долей миллиметра до 
нескольких миллиметров. Радиус перехода представляет собой случайную величину [2], может изменяться в 
достаточно широких пределах вдоль сварного соединения. Этот факт является причиной существенной 
нестабильности по длине шва значений концентрации напряжений, что необходимо учитывать при моделировании 
кинетики разрушения сварных соединений, работающих при переменных нагрузках. 

Для осуществления моделирования кинетики развития разрушения разработан алгоритм имитационного 
моделирования, который включает два этапа. 

На этапе 1 моделируется случайный процесс многоочагового зарождения трещин. С этой целью: 

1. Сварной шов разбивают на отдельные участки [; с одинаковыми значениями коэффициентов концентрации 
напряжений в пределах каждого участка в соответствии со статистическими характеристиками распределения 
значений длин участков. Положение участков вдоль сварного соединения Х; определяют с помощью датчика 
случайных чисел. 
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2. Для каждого участка в соответствии с кривой распределения номинальных напряжений (например, рис. 1) 
задают уровень номинальных напряжений о; (детерминированная величина). 

3. Используя полученную экспериментально функцию распределения радиусов, с помощью датчика 
случайных чисел для каждого Ё-того участка устанавливают значение радиуса и. 

4. Подставляя в выражение (2) значение радиуса, для каждого Г-того участка производят расчет значений 
теоретического (без учета развития пластических деформаций) коэффициента концентрации напряжений @. с учетом 
других необходимых параметров геометрии профиля сечения сварного соединения, которые можно принять 
неизменными вдоль всего шва (толщина стенки, размер катета и др.). 

5. Расчеты показывают, что на каждом участке поле главных напряжений имеет явно выраженный вектор, 
направленный перпендикулярно линии шва. Вторые два компонента имеют значительно меньшую величину по 
сравнению с первым, поэтому для упрощения в расчете напряженное состояние рассматривали как одноосное. 
Принимая во внимание, что металл в вершине концентратора находится в упругом состоянии, расчет коэффициента 
концентрации деформаций проводили как для одноосного растяжения с использованием закона Гука. Для каждого 
участка определяют амплитуду локальных деформаций на основе принятых значений номинальных напряжений (шаг 
2) и значений коэффициента концентрации (шаг 4). Т. к. локальные напряжения могут превышать предел текучести, 


то для расчета амплитуды деформаций используется уравнение Нейбера 


2 


2 2 
Ч = ар +, (3) 


где а, — коэффициент концентрации напряжений в упругой области; а, — коэффициент концентрации напряжений 
упруго-пластической области; а, — коэффициент концентрации деформаций. Уравнение (3) позволяет определить 
амплитуду локальных деформаций =а, необходимую для расчета длительности периода зарождения трещины. 

6. Для каждого Г-того участка рассчитывают число циклов до появления трещины в зависимости от 
амплитуды деформаций на {том участке. Для расчета используется модифицированное уравнение Басквина- 
Мэнсона-Коффина [6-8]: 





1 100 с'1/ Е" 
= 1 = , 4 
[8 п 1+и* пто0б— у с' | (1+») -- 
(4) + т 
1-* © (1+) 


6 
где Е" — модуль упругости; {® — относительное сужение; 0", — предел выносливости на базе 10°; оъ — предел 
прочности; т, — характеристика металла. Параметры Е", {", ; о принимаются в соответствии со стандартом или 
техническими условиями для конкретной стали. Значение показателя степени т, для малоуглеродистых 


низколегированных сталей с пределом прочности 0% = 300-700 МПа принимается равным 0,5. 

Расчеты, выполненные на этапе [ имитационного моделирования, позволяют создать базу Б; чисел циклов №; 
до зарождения трещин на каждом Ё-том участке вдоль сварного соединения. В качестве критерия зарождения 
разрушения принято появление полуэллиптической трещины длиной, равной длине участка [, глубиной В, = 0,5мм. 

На этапе П модель имитирует процесс развития многоочагового разрушения на основе уравнения Пэриса [12]: 

== С(Ак,)", (5) 
АМ 
где Г — скорость изменения глубины и длины полуэллиптической трещины в процессе дальнейшего циклического 
нагружения конструкции; АК; — размах коэффициента интенсивности напряжений; С и т — постоянные, 
соответствующие выбранному материалу. 

В отличие от роста одиночных трещин в данном случае слияние соседних трещин увеличивает длину фронта 
трещины на поверхности, что приводит к увеличению размаха коэффициента интенсивности напряжений и скорости 
роста трещины. 

С учетом случайного процесса накопления трещин вдоль сварного соединения, их подрастания и слияния 
фронтов соседних трещин, процесс имитационного моделирования развития разрушения включал в себя следующую 
последовательность действий: 

1. На основе базы данных числа циклов Б, до появления трещин определяется участок, на котором появились 
трещины при минимальном значении числа циклов №, и длина этого участка Г.,. Принимая для этого участка исходные 
размеры полуэллиптической трещины ‚= и Ь-0,5мм, производится расчет размаха коэффициента интенсивности 
напряжений. 

2. Значение числа циклов нагружения увеличивается на заданный шаг АМ и на основании уравнения Пэриса 
(5) рассчитываются новые размеры трещин, имевшихся на предыдущем шаге. 

3. Далее, используя базу Б;, определяют на каких участках сварного соединения должны появиться новые 
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Если фронты новых трещин не соприкасаются с существующими, то им назначают исходные размеры /Д, = / 
ир, = 0,5мм. 

Если фронты новых трещин соприкасаются с существующими фронтами, то осуществляют их слияние, при 
этом их длину на поверхности принимают равной сумме объединенной трещины, а в качестве глубины принимают 
наибольшую глубину из объединенных трещин. После этого производится сортировка базы Б‚, с присвоением новых 
номеров участков сварного соединения и размеров существующих на этом шаге трещин. Таким образом, получают 
откорректированную базу Б”, где и>1 — номер шага моделирования многоочагового развития разрушения. 

4. Используя базу Бг, увеличивают номер шага п+1 и повторяют операции, предусмотренные в пунктах 2—4 до 
тех пор, пока какая либо из трещин достигнет величины предельного размера. 

Предельный размер трещины устанавливают исходя из следующих соображений. Трещина предельного 
размера не должна нарушать эксплуатацию конструкции в штатном режиме, а визуальный контроль должен позволять 
выявить трещины при диагностике крупногабаритных конструкций. При обследовании парка мобильных кранов в 
новороссийском порту выяснилось, что минимальный размер обнаруживаемых при визуальном контроле трещин 
составляет 130 мм на поверхности. 

5. В связи с тем, что изменение геометрических параметров шва носит случайный характер, для получения 
статистически значимых результатов модель рассматривают многократно с измененным состоянием датчика 
случайных чисел. 

6. Расчет заканчивают статистической обработкой результатов моделирования с целью определения наиболее 
вероятного места появления трещин предельного размера и числа циклов нагружения. 

Выводы. Имитационное моделирование многоочагового разрушения позволяет: 

1. Прогнозировать периодичность проведения диагностики состояния конструкций; 

2. Упростить выявление мест повреждения; 

3. Обосновать технологические и конструктивные варианты ремонта и усиления рамной конструкции; 

4. Принимать решение о необходимости внесения конструктивных изменений для перераспределения 
силовых потоков в узле в случае обнаружения участков с локальным повышением концентрации напряжений. 
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